Origin and control of spin currents in a magnetic triplet Josephson junction 
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We study the appearance of a Josephson spin current in a model triplet superconductor junction 
with a magnetically-active tunnelling barrier. We find three distinct mechanisms for producing a 
spin current, and we provide a detailed discussion of the symmetry properties and the physical 
origins of each. By combining these three basic mechanisms, we find that it is possible to exercise 
fine control over the spin currents. In particular, we show that unlike the charge current, the spin 
currents on either side of the barrier need not be identical. 
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The fabrication and study of devices that combine 
triplet superconductivity and magnetism in a controlled 
manner will undoubtedly lead to new insights about the 
unique interplay of these two phases, and possibly also 
novel spintronic devices. Much attention has therefore 
been directed at developing the theoretical understand- 
ing of such systems [H, US S, S 0] • These efforts have al- 
ready revealed many unconventional effects in the trans- 
port properties of such devices: for example, the Joseph- 
son charge current through a triplet-superconductor- 
ferromagnet-triplet-superconductor (TFT) junction de- 
pends not only upon the magnitude (as in a singlet junc- 
tion, see Ref. 0) but also upon the orientation of the 
ferromagnetic moment at the barrier 0, 0] . The origin 
of this behaviour is the interaction of the barrier mo- 
ment with the spin-structure of the triplet state. De- 
spite the clear importance of the spin, only a few works 
have considered the appearance of a spontaneous spin 
current in materials where the pairing has a triplet com- 
ponent 0, H 1, i, . 

In this letter we study the appearance of a Josephson 
spin current in a minimal model of a TFT junction [H,[6j]. 
We find three general mechanisms for producing a spin 
current: spin-filtering, misalignment of the triplet vec- 
tor order parameters {d- vectors), and spin- flipping off a 
transverse barrier moment. The symmetry properties of 
the spin currents due to each mechanism are analyzed in 
detail, and the physical origins are understood by exam- 
ining the contributions to the charge and spin currents 
from each spin sector. Unlike the Josephson charge cur- 
rent, we show that the Josephson spin currents on each 
side of the barrier can be different. By judicious choice 
of the tunnelling barrier and the d- vector alignment, we 
can exercise fine control over the magnitude, direction 
and orientation in spin space of the spin currents on ei- 
ther side of the junction. 

The Hamiltonian describing the TFT junction is given 
by H — J TL(z, z')dzdz' with Hamiltonian density 
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FIG. 1: (color online) Schematic diagram of the model junc- 
tion studied in this work. 
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where v = L(R) for z',z < (z',z > 0). tpl{z) (il) a (z)) is 
the fermionic creation (annihilation) operator for a spin-<7 
particle with co-ordinate z. The quasiparticles have ef- 
fective mass m and chemical potential \i in the supercon- 
ductors. The superconducting gap is antisymmetric with 
respect to particle interchange, i.e. A(z,z r ) — ~A(z',z). 
As in Ref. d and Ref. @, we assume p 2 -pairing orbitals in 
both superconductors. The gap obeys a BCS tempera- 
ture dependence, and has T = magnitude Ao- The con- 
densates differ by a phase (p. The d- vector on the v side 
of the junction is parameterized as d v — (cos 9 V , sin 9 V , 0) 
where we take 9^ = and Or = 9 for definitcness (i.e. 
the left d- vector defines the x-axis). As the d- vectors 
have no z-component, this is the mutual spin orienta- 
tion for pairing in the two superconductors, i.e. they are 
both in an equal-spin-pairing state with respect to the 
z-axis. At the barrier we have potential scattering by 
Up(z) = UpS(z) and magnetic scattering by the moment 
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M(z,z') = (M± cos a, M± sin a, M\\)5(z)8{z'); <x Q/3 de- 
notes the vector of Pauli matrices. In what follows, we 
set h = 1 and quote values of M±, My and Up in units 
of m/kp where kp is the Fermi momentum. A schematic 
diagram of the junction is shown in Fig. |T]). 

We construct the 4x4 Nambu-spin space Green's func- 
tion Q (z , z' , f2; ito n ) of the junction in terms of quasi- 
classical scattering wavef unctions. A full account of this 
method for singlet superconductors is given in Ref. I 111 : 
the generalization to the triplet case is straight-forward. 
For simplicity, we assume spatially-constant order pa- 
rameters within the two superconductors. This should 
not qualitatively alter our results. Proceeding from the 
continuity equations, we define the Josephson charge (I") 
and //-component of the spin „) currents on the z^-side 
of the junction 
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where J n indicates an integration over half the spherical 
Fermi surface, f2 is the solid angle and Q u (z, z', O; icu n ) 
is the Green's function for z,z' > (<)0 when v = R(L). 
In the absence of chiral symmetry-breaking, the currents 
only flow along the z-axis. 

The most direct way to produce a spin current is by the 
spin- filtering effect [g]. For a spin-filter between equal- 
spin-pairing superconductors, we require both a finite Up 
and a component of the magnetic moment M in the 
spin quantization plane of each superconductor. The 
component of M in the spin quantization planes then 
defines a preferred direction, with the spin projections 
along this axis not mixed by the scattering at the bar- 
rier. Rather, we find spin-dependent effective scattering 
potentials Up ± \M\. The transparency of the barrier 
is therefore different in each spin sector, favouring the 
transmission of one spin species of Cooper pair over the 
other, i.e. the barrier selectively filters out one spin state, 
see Fig. ©(a). By time- reversal symmetry and mirror 
reflection of the junction in the x-z plane, we find that 
the spin current is antisymmetric in both <p and M. We 
now specialize to the case when M — Miiz, as this re- 
mains a spin-filter for arbitrary 9. The spin structure 
of the Josephson currents are best appreciated by isolat- 
ing the contributions to the total charge and spin cur- 
rents from each spin sector: we define the spin-er charge 
current as I" a = (I" — a2el v s z )/2, and the spin-cr spin 
current Ig Za = —al^ a /2e. For the spin-filter, we find 
that for Up, M\\ > the magnitude of 7 Cl f is greater than 
I c _l [Fig. ©(a)]. As such, Ig Z ^ and I" , , do not cancel, 
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FIG. 2: (color online) Cartoon representation of the three 
basic mechanisms for producing a spin current, (a) Spin- 
filter: spin-dependent effective barrier heights favour trans- 
mission through the spin-f sector over the spin- J. sector, (b) 
d- vector misalignment: each spin-sector sees a different effec- 
tive phase difference <j> — oQ. (c) Spin-flipping: electron-like 
quasiparticles spin-flipped by a transverse moment acquire a 
spin-dependent phase factor era. 



giving us the finite spin current shown in Fig. I[5])(b). 

Asano has shown that a spin current is also gener- 
ated in a non-magnetic junction by misaligning the d- 
vectors of the superconductors 0]. This produces a gra- 
dient of the order parameter in spin space, driving a 
spin current in analogy to the situation in superfluid 
3 He 12j. For the TFT junction studied here, this results 
in distinct phase differences in each of the spin chan- 
nels, 4> a = <j) — a9 0, [ijj], as illustrated schematically 
in Fig. (0)(b). We show the charge currents through each 
spin sector in Fig. ©(c): they are 29 out of phase with 
respect to one another but otherwise identical. Although 
each spin sector has a finite current flowing through it at 
= 0, these cancel and so the total charge current is van- 
ishing; this is not true of the spin current, demonstrating 
that the spin current can flow in the absence of a charge 
current [Fig. ©(d)]. Unlike the spin- filter, the spin cur- 
rent produced by the d- vector misalignment is symmetric 
with respect to </>, which follows from time-reversal sym- 
metry. When d- vector misalignment is combined with a 
spin-filter, the charge currents do not cancel at 6== 0, 7r 
and the junction is then in a fractional state @, [3|- 

Both the spin-filter and the tZ-vector misalignment 
mechanisms preserve the spin of the tunnelling Cooper 
pairs. We hence find that, like the charge current, the 
spin currents on either side of the junction are the same. 
This is not the case, however, for the third mechanism, 
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FIG. 3: (color online) Current vs phase relations for the three 
basic mechanisms, (a) Charge current and (b) ^-component of 
spin current for a spin-filter barrier (My = 0.5, Up = 1, other 
parameters vanishing), decomposed into contributions from 
each spin sector, (c) Charge current and (d) spin current for 
a non-magnetic barrier and with misaligned d-vectors of the 
two superconductors (Up = 1, 9 = 0.2-7T, other parameters 
vanishing), (e) Charge current and (f) spin current on the 
LHS of a barrier with a transverse moment (M± = 0.5, a — 
0.3-7T, other parameters vanishing). The spin current on the 
RHS is identical in magnitude but has opposite sign. In all 
figures we set T = 0.4T C . 



where the barrier magnetic moment has a component 
outside the quantization plane of the two superconduc- 
tors. We consider the least complicated such scenario: we 
have aligned d-vectors, and there is only a transverse mo- 
ment at the junction (i.e. M± ^ 0, M\\ = U P = 6 = 0). 
We find that only the z-component of the spin current 
is non- vanishing for < a < n/2. Although the magni- 
tude of the spin current is the same on each side of the 
junction, the sign is reversed, i.e. the v — L and v = R 
spin currents flow in opposite directions. This remark- 
able result is confirmed by a symmetry analysis: treating 
the spin current as a function only of the angle a and 
the phase <fi, we find by time-reversal and inversion sym- 
metry that I^ z {a,4>) = -I-^{a,<j>) [15]. For M ± = 0.5 
and a = 0.3ir, we show the charge current and its de- 



composition into spin-| and spin-! components on each 
side of the barrier in Fig. ©(e), while the z-component 
of the spin currents on the LHS are shown in Fig. ©(f). 



The physical interpretation of this behaviour follows by 
examining the charge currents through each spin sector: 
from Fig. ©(e), we observe that the current through the 
spin-T (!) sector on the LHS of the junction is equal to the 
current through the spin-! (!) sector on the RHS. That is, 
the spin of some of the tunnelling Cooper pairs is flipped 
by the transverse moment; this spin-relaxation process 
prevents a spin-accumulation occurring at the barrier. 
When a spin-cr electron-like (hole-like) quasiparticle tun- 
nelling between two normal regions undergoes a spin-flip 
at the tunnelling barrier, the transmitted and reflected 
waves acquire a phase shift of +(— )aa. Although the 
spin components of the charge current indicate that the 
situation is somewhat more complicated for tunnelling 
between two superconductors, this spin-dependent phase 
shift is ultimately the origin of the spin current. This 
conclusion is immediately apparent from the presence 
of the normal-normal scattering matrix in the bound- 
ary conditions for the cmasiclassical Green's function at 
a spin-active interface [l6( . We illustrate the spin-flipping 
mechanism schematically in Fig. (0(c). 

We have now discussed in detail the three basic mech- 
anisms for the production of a spin current in the TFT 
junction. We have already seen that the sign of the spin 
current need not be the same on either side of the junc- 
tion. By combining the three mechanisms, however, we 
demonstrate that we can also vary the magnitude and 
spin-orientation of the current. 

Upon combining the transverse moment with a spin- 
filter (Mx,M«,Up ^ 0, 9 = 0), we find two significant 
changes: a y-component of the spin current appears in 
addition to the z-component, and the magnitude of the 
spin currents on either side of the barrier are no longer 
the same for < a < ir/2. We show the spin currents 
on the LHS and RHS of such a barrier in Fig. ©(a) 
and (b) respectively. We consider first the origins of the 
y-component. For a =/= 0,7r, this follows immediately 
from the spin-filter effect: for aligned d-vectors, the y- 
axis lies within the quantization plane of both supercon- 
ductors. As the magnetization lying in the y-z plane has 
a finite component along the y-axis, the resulting spin- 
filter selects a direction with a finite y-component, and 
so I" y ^ 0. This is the only magnetic interaction at the 
barrier when a = n/2, and the spin current is then the 
same on either side of the junction. For < a < n/2, 
however, there is a finite component of the magnetiza- 
tion along each axis: the scattering off the components 
transverse to the y-axis and off those transverse to the 
z-axis also produces a y- and a z-component of the spin 
current respectively. As above, the spin currents pro- 
duced by the spin-flipping flow in opposite directions on 
either side of the junction. The total spin current on 
each side of the junction is then the sum of a contri- 
bution which is symmetric and a contribution which is 
antisymmetric with respect to v. As the former is anti- 
symmetric in <p while the latter is symmetric, and both 
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FIG. 4: (color online) /i-component of the spin currents from 
combinations of the three mechanisms. Spin currents on the 
(a) LHS and (b) RHS of a junction combining spin-filtering 
with spin-flipping (Up = 1, Mm = M± = 0.5, a — 0.3ir, other 
parameters vanishing). Spin currents on the (c) LHS and (d) 
RHS of a junction combining a transverse barrier moment and 
misaligned d-vectors (M± = 0.5, a = 0.5-7T, 9 — 0.2ir, other 
parameters vanishing). In all figures we set T = 0.4T C . 



tween the v — L and v = R z-component of the spin 
current: this is strongly enhanced at = 7r on the RHS 
over that on the LHS, while its value at 4> = is much 
the same in the two cases. 

In conclusion, we have used the quasiclassical method 
to study the interaction of triplet superconductivity and 
magnetism in a model Josephson junction. Our main re- 
sult is the classification of the three distinct mechanisms 
for the production of a spin current. The signature prop- 
erties of each of these mechanisms are given intuitive 
physical explanations by decomposing the total charge 
and spin currents into contributions from each spin sec- 
tor [Fig. ©]. We demonstrate that in the presence of 
a spin-relaxing moment at the barrier, neither the mag- 
nitude, direction or spin orientation of the spin currents 
on each side of the junction need be identical. Our work 
theoretically opens the route for the inclusion of triplet 
Josephson junctions in spintronic applications. 

The authors thank C. Iniotakis, Y. Asano, K. Benne- 
mann, B. Kastening, D.K. Morr and F.S. Nogueira for 
useful discussions. The authors also thank J. Sirker for 
his careful reading of the manuscript. P.M.R.B. grate- 
fully acknowledges the hospitality of A. Simon. 



are 2-7r-periodic in the phase difference, we have the re- 
lation I v s z(y) ((j)) = -I~z( y )(2K - <P), which can clearly be 
seen in Fig. Q. 

We finally discuss the combination of a spin non- 
conserving moment with a misalignment of the d- vectors. 
For a finite 6 ^ 0, ir, the quantization planes of the two 
superconductors only intersect along the z-axis: any bar- 
rier moment which has a component transverse to the 
z-axis is not spin-conserving for the tunnelling Cooper 
pairs. We therefore consider the case where, in the ab- 
sence of the d- vector misalignment, the magnetic scatter- 
ing produces a y-component of the spin current. When 
the d- vectors are misaligned, a y-component of the spin 
current still flows on the LHS; on the RHS, however, 
we have an additional x-component of current, which is 
related to the y-component by If = twi(6)lf y . By us- 
ing a magnetic moment at the barrier, we are therefore 
able to produce a spin current on the RHS with a com- 
ponent along the axis perpendicular to the quantization 
plane in the LHS superconductor. This situation is dis- 
played in Fig. (H])(c) and (d), which respectively show 
the spin currents on the LHS and RHS of a spin-filter for 
the y-component of spin (Mj_, Up ^ 0, a = ir/2). Note 
that the finite z-component of the spin current is com- 
posed of contributions due to the misaligned ti-vectors 
and also from the spin-flip scattering. Unlike the previ- 
ous cases, there is no simple relationship between the spin 
currents on either side of the barrier, which is reflected 
in Fig. (d])(c) and (d) by the considerable differences be- 
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